Posttranscriptional regulation of FAN1 by miR-124-3p at rs3512 underlies onset-delaying genetic modification in Huntington's disease.

Many Mendelian disorders, such as Huntington's disease (HD) and spinocerebellar ataxias, arise from expansions of CAG trinucleotide repeats. Despite the clear genetic causes, additional genetic factors may influence the rate of those monogenic disorders. Notably, genome-wide association studies discovered somewhat expected modifiers, particularly mismatch repair genes involved in the CAG repeat instability, impacting age at onset of HD. Strikingly, FAN1, previously unrelated to repeat instability, produced the strongest HD modification signals. Diverse FAN1 haplotypes independently modify HD, with rare genetic variants diminishing DNA binding or nuclease activity of the FAN1 protein, hastening HD onset. However, the mechanism behind the frequent and the most significant onset-delaying FAN1 haplotype lacking missense variations has remained elusive. Here, we illustrated that a microRNA acting on 3'-UTR (untranslated region) SNP rs3512, rather than transcriptional regulation, is responsible for the significant FAN1 expression quantitative trait loci signal and allelic imbalance in FAN1 messenger ribonucleic acid (mRNA), accounting for the most significant and frequent onset-delaying modifier haplotype in HD. Specifically, miR-124-3p selectively targets the reference allele at rs3512, diminishing the stability of FAN1 mRNA harboring that allele and consequently reducing its levels. Subsequent validation analyses, including the use of antagomir and 3'-UTR reporter vectors with swapped alleles, confirmed the specificity of miR-124-3p at rs3512. Together, these findings indicate that the alternative allele at rs3512 renders the FAN1 mRNA less susceptible to miR-124-3p-mediated posttranscriptional regulation, resulting in increased FAN1 levels and a subsequent delay in HD onset by mitigating CAG repeat instability.

Somatic CAG repeat instability in intermediate alleles of the HTT gene and its potential association with a clinical phenotype.

Huntington disease (HD) is a neurodegenerative disorder caused by ≥36 CAGs in the HTT gene. Intermediate alleles (IAs) (27-35 CAGs) are not considered HD-causing, but their potential association with neurocognitive symptoms remains controversial. As HTT somatic CAG expansion influences HD onset, we hypothesised that IAs are somatically unstable, and that somatic CAG expansion may drive phenotypic presentation in some IA carriers. We quantified HTT somatic CAG expansions by MiSeq sequencing in the blood DNA of 164 HD subjects and 191 IA (symptomatic and control) carriers, and in the brain DNA of a symptomatic 33 CAG carrier. We also performed genotype-phenotype analysis. The phenotype of symptomatic IA carriers was characterised by motor (85%), cognitive (27%) and/or behavioural (29%) signs, with a late (58.7 ± 18.6 years), but not CAG-dependent, age at onset. IAs displayed somatic expansion that were CAG and age-dependent in blood DNA, with 0.4% and 0.01% of DNA molecules expanding by CAG and year, respectively. Somatic expansions of +1 and +2 CAGs were detected in the brain of the individual with 33 CAGs, with the highest expansion frequency in the putamen (10.3%) and the lowest in the cerebellum (4.8%). Somatic expansion in blood DNA was not different in symptomatic vs. control IA carriers. In conclusion, we show that HTT IAs are somatically unstable, but we found no association with HD-like phenotypes. It is plausible, however, that some IAs, close to the HD pathological threshold and with a predisposing genetic background, could manifest with neurocognitive symptoms.

Modification of Huntington's disease by short tandem repeats.

Expansions of glutamine-coding CAG trinucleotide repeats cause a number of neurodegenerative diseases, including Huntington's disease and several of spinocerebellar ataxias. In general, age-at-onset of the polyglutamine diseases is inversely correlated with the size of the respective inherited expanded CAG repeat. Expanded CAG repeats are also somatically unstable in certain tissues, and age-at-onset of Huntington's disease corrected for individual HTT CAG repeat length (i.e. residual age-at-onset), is modified by repeat instability-related DNA maintenance/repair genes as demonstrated by recent genome-wide association studies. Modification of one polyglutamine disease (e.g. Huntington's disease) by the repeat length of another (e.g. ATXN3, CAG expansions in which cause spinocerebellar ataxia 3) has also been hypothesized. Consequently, we determined whether age-at-onset in Huntington's disease is modified by the CAG repeats of other polyglutamine disease genes. We found that the CAG measured repeat sizes of other polyglutamine disease genes that were polymorphic in Huntington's disease participants but did not influence Huntington's disease age-at-onset. Additional analysis focusing specifically on ATXN3 in a larger sample set (n = 1388) confirmed the lack of association between Huntington's disease residual age-at-onset and ATXN3 CAG repeat length. Additionally, neither our Huntington's disease onset modifier genome-wide association studies single nucleotide polymorphism data nor imputed short tandem repeat data supported the involvement of other polyglutamine disease genes in modifying Huntington's disease. By contrast, our genome-wide association studies based on imputed short tandem repeats revealed significant modification signals for other genomic regions. Together, our short tandem repeat genome-wide association studies show that modification of Huntington's disease is associated with short tandem repeats that do not involve other polyglutamine disease-causing genes, refining the landscape of Huntington's disease modification and highlighting the importance of rigorous data analysis, especially in genetic studies testing candidate modifiers.

Spanish HTT gene study reveals haplotype and allelic diversity with possible implications for germline expansion dynamics in Huntington disease.

We aimed to determine the genetic diversity and molecular characteristics of the Huntington disease (HD) gene (HTT) in Spain. We performed an extended haplotype and exon one deep sequencing analysis of the HTT gene in a nationwide cohort of population-based controls (n = 520) and families with symptomatic individuals referred for HD genetic testing. This group included 331 HD cases and 140 carriers of intermediate alleles. Clinical and family history data were obtained when available. Spanish normal alleles are enriched in C haplotypes (40.1%), whereas A1 (39.8%) and A2 (31.6%) prevail among intermediate and expanded alleles, respectively. Alleles ≥ 50 CAG repeats are primarily associated with haplotypes A2 (38.9%) and C (32%), which are also present in 50% and 21.4%, respectively, of HD families with large intergenerational expansions. Non-canonical variants of exon one sequence are less frequent, but much more diverse, in alleles of ≥27 CAG repeats. The deletion of CAACAG, one of the six rare variants not observed among smaller normal alleles, is associated with haplotype C and appears to correlate with larger intergenerational expansions and early onset of symptoms. Spanish HD haplotypes are characterized by a high genetic diversity, potentially admixed with other non-Caucasian populations, with a higher representation of A2 and C haplotypes than most European populations. Differences in haplotype distributions across the CAG length range support differential germline expansion dynamics, with A2 and C showing the largest intergenerational expansions. This haplotype-dependent germline instability may be driven by specific cis-elements, such as the CAACAG deletion.

Repeat Detector: versatile sizing of expanded tandem repeats and identification of interrupted alleles from targeted DNA sequencing.

Targeted DNA sequencing approaches will improve how the size of short tandem repeats is measured for diagnostic tests and preclinical studies. The expansion of these sequences causes dozens of disorders, with longer tracts generally leading to a more severe disease. Interrupted alleles are sometimes present within repeats and can alter disease manifestation. Determining repeat size mosaicism and identifying interruptions in targeted sequencing datasets remains a major challenge. This is in part because standard alignment tools are ill-suited for repetitive and unstable sequences. To address this, we have developed Repeat Detector (RD), a deterministic profile weighting algorithm for counting repeats in targeted sequencing data. We tested RD using blood-derived DNA samples from Huntington's disease and Fuchs endothelial corneal dystrophy patients sequenced using either Illumina MiSeq or Pacific Biosciences single-molecule, real-time sequencing platforms. RD was highly accurate in determining repeat sizes of 609 blood-derived samples from Huntington's disease individuals and did not require prior knowledge of the flanking sequences. Furthermore, RD can be used to identify alleles with interruptions and provide a measure of repeat instability within an individual. RD is therefore highly versatile and may find applications in the diagnosis of expanded repeat disorders and in the development of novel therapies.

A probable cis-acting genetic modifier of Huntington disease frequent in individuals with African ancestry.

Huntington disease (HD)is a dominantly inherited neurodegenerative disorder caused by the expansion of a polyglutamine encoding CAG repeat in the huntingtin gene. Recently, it has been established that disease severity in HD is best predicted by the number of pure CAG repeats rather than total glutamines encoded. Along with uncovering DNA repair gene variants as trans-acting modifiers of HD severity, these data reveal somatic expansion of the CAG repeat as a key driver of HD onset. Using high-throughput DNA sequencing, we have determined the precise sequence and somatic expansion profiles of the HTT repeat tract of 68 HD-affected and 158 HD-unaffected African ancestry individuals. A high level of HTT repeat sequence diversity was observed, with three likely African-specific alleles identified. In the most common disease allele (30 out of 68), the typical proline-encoding CCGCCA sequence was absent. This CCGCCA-loss disease allele was associated with an earlier age of diagnosis of approximately 7.1 years and occurred exclusively on haplotype B2. Although somatic expansion was associated with an earlier age of diagnosis in the study overall, the CCGCCA-loss disease allele displayed reduced somatic expansion relative to the typical HTT expansions in blood DNA. We propose that the CCGCCA loss occurring on haplotype B2 is an African cis-acting modifier that appears to alter disease diagnosis of HD through a mechanism that is not driven by somatic expansion. The assessment of a group of individuals from an understudied population has highlighted population-specific differences that emphasize the importance of studying genetically diverse populations in the context of disease.

Exome sequencing of individuals with Huntington's disease implicates FAN1 nuclease activity in slowing CAG expansion and disease onset.

The age at onset of motor symptoms in Huntington's disease (HD) is driven by HTT CAG repeat length but modified by other genes. In this study, we used exome sequencing of 683 patients with HD with extremes of onset or phenotype relative to CAG length to identify rare variants associated with clinical effect. We discovered damaging coding variants in candidate modifier genes identified in previous genome-wide association studies associated with altered HD onset or severity. Variants in FAN1 clustered in its DNA-binding and nuclease domains and were associated predominantly with earlier-onset HD. Nuclease activities of purified variants in vitro correlated with residual age at motor onset of HD. Mutating endogenous FAN1 to a nuclease-inactive form in an induced pluripotent stem cell model of HD led to rates of CAG expansion similar to those observed with complete FAN1 knockout. Together, these data implicate FAN1 nuclease activity in slowing somatic repeat expansion and hence onset of HD.

Genetic modifiers of Huntington disease differentially influence motor and cognitive domains.

Genome-wide association studies (GWASs) of Huntington disease (HD) have identified six DNA maintenance gene loci (among others) as modifiers and implicated a two step-mechanism of pathogenesis: somatic instability of the causative HTT CAG repeat with subsequent triggering of neuronal damage. The largest studies have been limited to HD individuals with a rater-estimated age at motor onset. To capitalize on the wealth of phenotypic data in several large HD natural history studies, we have performed algorithmic prediction by using common motor and cognitive measures to predict age at other disease landmarks as additional phenotypes for GWASs. Combined with imputation with the Trans-Omics for Precision Medicine reference panel, predictions using integrated measures provided objective landmark phenotypes with greater power to detect most modifier loci. Importantly, substantial differences in the relative modifier signal across loci, highlighted by comparing common modifiers at MSH3 and FAN1, revealed that individual modifier effects can act preferentially in the motor or cognitive domains. Individual components of the DNA maintenance modifier mechanisms may therefore act differentially on the neuronal circuits underlying the corresponding clinical measures. In addition, we identified additional modifier effects at the PMS1 and PMS2 loci and implicated a potential second locus on chromosome 7. These findings indicate that broadened discovery and characterization of HD genetic modifiers based on additional quantitative or qualitative phenotypes offers not only the promise of in-human validated therapeutic targets but also a route to dissecting the mechanisms and cell types involved in both the somatic instability and toxicity components of HD pathogenesis.

Blood meal analysis of tsetse flies ( Glossina pallidipes: Glossinidae) reveals higher host fidelity on wild compared with domestic hosts.

Background: Changes in climate and land use can alter risk of transmission of parasites between domestic hosts and wildlife, particularly when mediated by vectors that can travel between populations. Here we focused on tsetse flies (genus Glossina), the cyclical vectors for both Human African Trypanosomiasis (HAT) and Animal African Trypanosomiasis (AAT). The aims of this study were to investigate three issues related to G. palldipes from Kenya: 1) the diversity of vertebrate hosts that flies fed on; 2) whether host feeding patterns varied in relation to type of hosts, tsetse feeding behaviour, site or tsetse age and sex; and 3) if there was a relationship between trypanosome detection and host feeding behaviours or host types. Methods: Sources of blood meals of Glossina pallidipes were identified by sequencing of the mitochondrial cytochrome b gene and analyzed in relationship with previously determined trypanosome detection in the same flies. Results: In an area dominated by wildlife but with seasonal presence of livestock (Nguruman), 98% of tsetse fed on single wild host species, whereas in an area including a mixture of resident domesticated animals, humans and wildlife (Shimba Hills), 52% of flies fed on more than one host species. Multiple Correspondence Analysis revealed strong correlations between feeding pattern, host type and site but these were resolved along a different dimension than trypanosome status, sex and age of the flies. Conclusions: Our results suggest that individual G. pallidipes in interface areas may show higher feeding success on wild hosts when available but often feed on both wild and domesticated hosts. This illustrates the importance of G. pallidipes as a vector connecting the sylvatic and domestic cycles of African trypanosomes.

Association Analysis of Chromosome X to Identify Genetic Modifiers of Huntington's Disease.

BACKGROUND: Huntington's disease (HD) is caused by an expanded (>35) CAG trinucleotide repeat in huntingtin (HTT). Age-at-onset of motor symptoms is inversely correlated with the size of the inherited CAG repeat, which expands further in brain regions due to somatic repeat instability. Our recent genetic investigation focusing on autosomal SNPs revealed that age-at-onset is also influenced by genetic variation at many loci, the majority of which encode genes involved in DNA maintenance/repair processes and repeat instability. OBJECTIVE: We performed a complementary association analysis to determine whether variants in the X chromosome modify HD. METHODS: We imputed SNPs on chromosome X for ∼9,000 HD subjects of European ancestry and performed an X chromosome-wide association study (XWAS) to test for association with age-at-onset corrected for inherited CAG repeat length. RESULTS: In a mixed effects model XWAS analysis of all subjects (males and females), assuming random X-inactivation in females, no genome-wide significant onset modification signal was found. However, suggestive significant association signals were detected at Xq12 (top SNP, rs59098970; p-value, 1.4E-6), near moesin (MSN), in a region devoid of DNA maintenance genes. Additional suggestive signals not involving DNA repair genes were observed in male- and female-only analyses at other locations. CONCLUSION: Although not genome-wide significant, potentially due to small effect size compared to the power of the current study, our data leave open the possibility of modification of HD by a non-DNA repair process. Our XWAS results are publicly available at the updated GEM EURO 9K website hosted at https://www.hdinhd.org/ for browsing, pathway analysis, and data download.

Approaches to Sequence the HTT CAG Repeat Expansion and Quantify Repeat Length Variation.

BACKGROUND: Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder caused by the expansion of the HTT CAG repeat. Affected individuals inherit ≥36 repeats and longer alleles cause earlier onset, greater disease severity and faster disease progression. The HTT CAG repeat is genetically unstable in the soma in a process that preferentially generates somatic expansions, the proportion of which is associated with disease onset, severity and progression. Somatic mosaicism of the HTT CAG repeat has traditionally been assessed by semi-quantitative PCR-electrophoresis approaches that have limitations (e.g., no information about sequence variants). Genotyping-by-sequencing could allow for some of these limitations to be overcome. OBJECTIVE: To investigate the utility of PCR sequencing to genotype large (>50 CAGs) HD alleles and to quantify the associated somatic mosaicism. METHODS: We have applied MiSeq and PacBio sequencing to PCR products of the HTT CAG repeat in transgenic R6/2 mice carrying ∼55, ∼110, ∼255 and ∼470 CAGs. For each of these alleles, we compared the repeat length distributions generated for different tissues at two ages. RESULTS: We were able to sequence the CAG repeat full length in all samples. However, the repeat length distributions for samples with ∼470 CAGs were biased towards shorter repeat lengths. CONCLUSION: PCR sequencing can be used to sequence all the HD alleles considered, but this approach cannot be used to estimate modal allele size or quantify somatic expansions for alleles ⪢250 CAGs. We review the limitations of PCR sequencing and alternative approaches that may allow the quantification of somatic contractions and very large somatic expansions.

DNA methylation study of Huntington's disease and motor progression in patients and in animal models.

Although Huntington's disease (HD) is a well studied Mendelian genetic disorder, less is known about its associated epigenetic changes. Here, we characterize DNA methylation levels in six different tissues from 3 species: a mouse huntingtin (Htt) gene knock-in model, a transgenic HTT sheep model, and humans. Our epigenome-wide association study (EWAS) of human blood reveals that HD mutation status is significantly (p < 10-7) associated with 33 CpG sites, including the HTT gene (p = 6.5 × 10-26). These Htt/HTT associations were replicated in the Q175 Htt knock-in mouse model (p = 6.0 × 10-8) and in the transgenic sheep model (p = 2.4 × 10-88). We define a measure of HD motor score progression among manifest HD cases based on multiple clinical assessments. EWAS of motor progression in manifest HD cases exhibits significant (p < 10-7) associations with methylation levels at three loci: near PEX14 (p = 9.3 × 10-9), GRIK4 (p = 3.0 × 10-8), and COX4I2 (p = 6.5 × 10-8). We conclude that HD is accompanied by profound changes of DNA methylation levels in three mammalian species.

Genetic Risk Underlying Psychiatric and Cognitive Symptoms in Huntington's Disease.

BACKGROUND: Huntington's disease (HD) is an inherited neurodegenerative disorder caused by an expanded CAG repeat in the HTT gene. It is diagnosed following a standardized examination of motor control and often presents with cognitive decline and psychiatric symptoms. Recent studies have detected genetic loci modifying the age at onset of motor symptoms in HD, but genetic factors influencing cognitive and psychiatric presentations are unknown. METHODS: We tested the hypothesis that psychiatric and cognitive symptoms in HD are influenced by the same common genetic variation as in the general population by 1) constructing polygenic risk scores from large genome-wide association studies of psychiatric and neurodegenerative disorders and of intelligence and 2) testing for correlation with the presence of psychiatric and cognitive symptoms in a large sample (n = 5160) of patients with HD. RESULTS: Polygenic risk score for major depression was associated specifically with increased risk of depression in HD, as was schizophrenia risk score with psychosis and irritability. Cognitive impairment and apathy were associated with reduced polygenic risk score for intelligence. CONCLUSIONS: Polygenic risk scores for psychiatric disorders, particularly depression and schizophrenia, are associated with increased risk of the corresponding psychiatric symptoms in HD, suggesting a common genetic liability. However, the genetic liability to cognitive impairment and apathy appears to be distinct from other psychiatric symptoms in HD. No associations were observed between HD symptoms and risk scores for other neurodegenerative disorders. These data provide a rationale for treatments effective in depression and schizophrenia to be used to treat depression and psychotic symptoms in HD.

A genetic association study of glutamine-encoding DNA sequence structures, somatic CAG expansion, and DNA repair gene variants, with Huntington disease clinical outcomes.

BACKGROUND: Huntington disease (HD) is caused by an unstable CAG/CAA repeat expansion encoding a toxic polyglutamine tract. Here, we tested the hypotheses that HD outcomes are impacted by somatic expansion of, and polymorphisms within, the HTT CAG/CAA glutamine-encoding repeat, and DNA repair genes. METHODS: The sequence of the glutamine-encoding repeat and the proportion of somatic CAG expansions in blood DNA from participants inheriting 40 to 50 CAG repeats within the TRACK-HD and Enroll-HD cohorts were determined using high-throughput ultra-deep-sequencing. Candidate gene polymorphisms were genotyped using kompetitive allele-specific PCR (KASP). Genotypic associations were assessed using time-to-event and regression analyses. FINDINGS: Using data from 203 TRACK-HD and 531 Enroll-HD participants, we show that individuals with higher blood DNA somatic CAG repeat expansion scores have worse HD outcomes: a one-unit increase in somatic expansion score was associated with a Cox hazard ratio for motor onset of 3·05 (95% CI = 1·94 to 4·80, p = 1·3 × 10-6). We also show that individual-specific somatic expansion scores are associated with variants in FAN1 (pFDR = 4·8 × 10-6), MLH3 (pFDR = 8·0 × 10-4), MLH1 (pFDR = 0·004) and MSH3 (pFDR = 0·009). We also show that HD outcomes are best predicted by the number of pure CAGs rather than total encoded-glutamines. INTERPRETATION: These data establish pure CAG length, rather than encoded-glutamine, as the key inherited determinant of downstream pathophysiology. These findings have implications for HD diagnostics, and support somatic expansion as a mechanistic link for genetic modifiers of clinical outcomes, a driver of disease, and potential therapeutic target in HD and related repeat expansion disorders. FUNDING: CHDI Foundation.

MSH3 modifies somatic instability and disease severity in Huntington's and myotonic dystrophy type 1.

The mismatch repair gene MSH3 has been implicated as a genetic modifier of the CAG·CTG repeat expansion disorders Huntington's disease and myotonic dystrophy type 1. A recent Huntington's disease genome-wide association study found rs557874766, an imputed single nucleotide polymorphism located within a polymorphic 9 bp tandem repeat in MSH3/DHFR, as the variant most significantly associated with progression in Huntington's disease. Using Illumina sequencing in Huntington's disease and myotonic dystrophy type 1 subjects, we show that rs557874766 is an alignment artefact, the minor allele for which corresponds to a three-repeat allele in MSH3 exon 1 that is associated with a reduced rate of somatic CAG·CTG expansion (P = 0.004) and delayed disease onset (P = 0.003) in both Huntington's disease and myotonic dystrophy type 1, and slower progression (P = 3.86 × 10-7) in Huntington's disease. RNA-Seq of whole blood in the Huntington's disease subjects found that repeat variants are associated with MSH3 and DHFR expression. A transcriptome-wide association study in the Huntington's disease cohort found increased MSH3 and DHFR expression are associated with disease progression. These results suggest that variation in the MSH3 exon 1 repeat region influences somatic expansion and disease phenotype in Huntington's disease and myotonic dystrophy type 1, and suggests a common DNA repair mechanism operates in both repeat expansion diseases.

Sodalis glossinidius presence in wild tsetse is only associated with presence of trypanosomes in complex interactions with other tsetse-specific factors.

BACKGROUND: Susceptibility of tsetse flies (Glossina spp.) to trypanosomes of both humans and animals has been associated with the presence of the endosymbiont Sodalis glossinidius. However, intrinsic biological characteristics of the flies and environmental factors can influence the presence of both S. glossinidius and the parasites. It thus remains unclear whether it is the S. glossinidius or other attributes of the flies that explains the apparent association. The objective of this study was to test whether the presence of Trypanosoma vivax, T. congolense and T. brucei are related to the presence of S. glossinidius in tsetse flies when other factors are accounted for: geographic location, species of Glossina, sex or age of the host flies. RESULTS: Flies (n = 1090) were trapped from four sites in the Shimba Hills and Nguruman regions in Kenya. Sex and species of tsetse (G. austeni, G. brevipalpis, G. longipennis and G. pallidipes) were determined based on external morphological characters and age was estimated by a wing fray score method. The presence of trypanosomes and S. glossinidius was detected using PCR targeting the internal transcribed spacer region 1 and the haemolysin gene, respectively. Sequencing was used to confirm species identification. Generalised Linear Models (GLMs) and Multiple Correspondence Analysis (MCA) were applied to investigate multivariable associations. The overall prevalence of trypanosomes was 42.1%, but GLMs revealed complex patterns of associations: the presence of S. glossinidius was associated with trypanosome presence but only in interactions with other factors and only in some species of trypanosomes. The strongest association was found for T. congolense, and no association was found for T. vivax. The MCA also suggested only a weak association between the presence of trypanosomes and S. glossinidius. Trypanosome-positive status showed strong associations with sex and age while S. glossinidius-positive status showed a strong association with geographic location and species of fly. CONCLUSIONS: We suggest that previous conclusions about the presence of endosymbionts increasing probability of trypanosome presence in tsetse flies may have been confounded by other factors, such as community composition of the tsetse flies and the specific trypanosomes found in different regions.

Molecular identification of different trypanosome species and subspecies in tsetse flies of northern Nigeria.

BACKGROUND: Animal African Trypanosomiasis (AAT) is caused by several species of trypanosomes including Trypanosoma congolense, T. vivax, T. godfreyi, T. simiae and T. brucei. Two of the subspecies of T. brucei also cause Human African Trypanosomiasis. Although some of them can be mechanically transmitted by biting flies; these trypanosomes are all transmitted by tsetse flies which are the cyclical vectors of Trypanosoma congolense, T. godfreyi, T. simiae and T. brucei. We present here the first report assessing the prevalence of trypanosomes in tsetse flies in Nigeria using molecular tools. METHODS: 488 tsetse flies of three species, Glossina palpalis palpalis, G. tachinoides and G. morsitans submorsitans were collected from Wuya, Niger State and Yankari National Park, Bauchi State in 2012. Trypanosomes were detected and identified using an ITS1 PCR assay on DNA purified from the 'head plus proboscis' (H + P) and abdomen (ABD) parts of each fly. RESULTS: T. vivax and T. congolense Savannah were the major parasites detected. Trypanosomes prevalence was 7.1 % in G. p. palpalis, 11.9 % in G. tachinoides and 13.5 % in G. m. submorsitans. Prevalences of T. congolense Savannah ranged from 2.5 to 6.7 % and of T. vivax were approximately 4.5 %. Trypanosoma congolense Forest, T. godfreyi and T. simiae were also detected in the site of Yankari. The main biological and ecological determinants of trypanosome prevalence were the fly sex, with more trypanosomes found in females than males, and the site, with T. congolense subspp. being more abundant in Yankari than in Wuya. As expected, the trypanosome species diversity was higher in Yankari National Park than in the more agricultural site of Wuya where vertebrate host species diversity is lower. CONCLUSIONS: Our results show that T. congolense Savannah and T. vivax are the main species of parasite potentially causing AAT in the two study sites and that Yankari National Park is a potential reservoir of trypanosomes both in terms of parasite abundance and species diversity.

Laboratory colonisation and genetic bottlenecks in the tsetse fly Glossina pallidipes.

BACKGROUND: The IAEA colony is the only one available for mass rearing of Glossina pallidipes, a vector of human and animal African trypanosomiasis in eastern Africa. This colony is the source for Sterile Insect Technique (SIT) programs in East Africa. The source population of this colony is unclear and its genetic diversity has not previously been evaluated and compared to field populations. METHODOLOGY/PRINCIPAL FINDINGS: We examined the genetic variation within and between the IAEA colony and its potential source populations in north Zimbabwe and the Kenya/Uganda border at 9 microsatellites loci to retrace the demographic history of the IAEA colony. We performed classical population genetics analyses and also combined historical and genetic data in a quantitative analysis using Approximate Bayesian Computation (ABC). There is no evidence of introgression from the north Zimbabwean population into the IAEA colony. Moreover, the ABC analyses revealed that the foundation and establishment of the colony was associated with a genetic bottleneck that has resulted in a loss of 35.7% of alleles and 54% of expected heterozygosity compared to its source population. Also, we show that tsetse control carried out in the 1990's is likely reduced the effective population size of the Kenya/Uganda border population. CONCLUSIONS/SIGNIFICANCE: All the analyses indicate that the area of origin of the IAEA colony is the Kenya/Uganda border and that a genetic bottleneck was associated with the foundation and establishment of the colony. Genetic diversity associated with traits that are important for SIT may potentially have been lost during this genetic bottleneck which could lead to a suboptimal competitiveness of the colony males in the field. The genetic diversity of the colony is lower than that of field populations and so, studies using colony flies should be interpreted with caution when drawing general conclusions about G. pallidipes biology.

Secondary contact and admixture between independently invading populations of the western corn rootworm, Diabrotica virgifera virgifera in Europe.

The western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is one of the most destructive pests of corn in North America and is currently invading Europe. The two major invasive outbreaks of rootworm in Europe have occurred, in North-West Italy and in Central and South-Eastern Europe. These two outbreaks originated from independent introductions from North America. Secondary contact probably occurred in North Italy between these two outbreaks, in 2008. We used 13 microsatellite markers to conduct a population genetics study, to demonstrate that this geographic contact resulted in a zone of admixture in the Italian region of Veneto. We show that i) genetic variation is greater in the contact zone than in the parental outbreaks; ii) several signs of admixture were detected in some Venetian samples, in a bayesian analysis of the population structure and in an approximate bayesian computation analysis of historical scenarios and, finally, iii) allelic frequency clines were observed at microsatellite loci. The contact between the invasive outbreaks in North-West Italy and Central and South-Eastern Europe resulted in a zone of admixture, with particular characteristics. The evolutionary implications of the existence of a zone of admixture in Northern Italy and their possible impact on the invasion success of the western corn rootworm are discussed.

Biological invasions in agricultural settings: insights from evolutionary biology and population genetics.

Invasion biology and agriculture are intimately related for several reasons and in particular because many agricultural pest species are recent invaders. In this article we suggest that the reconstruction of invasion routes with population genetics-based methods can address fundamental questions in ecology and practical aspects of the management of biological invasions in agricultural settings. We provide a brief description of the methods used to reconstruct invasion routes and describe their main characteristics. In particular, we focus on a scenario--the bridgehead invasion scenario --which had been overlooked until recently. We show that this scenario, in which an invasive population is the source of other invasive populations, is evolutionarily parsimonious and may have played a crucial role in shaping the distribution of many recent agricultural pests.